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PREFACE 
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1.0  INTRODUCTION 

Interest  in  highly  maneuverable,  aerodynamically  controlled  missiles  has  increased  in 
recent  years.  Evaluation  of  mission  capability  for  missiles  of  this  category  relies  primarily  on 
the  missile  aerodynamics.  It  is  well  known  that  the  missile  static  aerodynamics  are 
predominant  in  a  mission  analysis;  what  is  not  known  is  the  importance  of  the  dynamic 
derivatives  in  such  an  analysis.  It  is  necessary  to  define  the  effects  of  these  parameters  so  that 
in  the  future,  an  experimental  test  apparatus  may  be  designed  for  their  measurement  and  so 
that  accurate  simulations  can  be  performed. 

This  subject  analysis  documents  the  effects  of  only  those  dynamic  derivatives  that  may 
significantly  affect  the  motion  of  highly  maneuverable  missiles.  The  derivatives  investigated 
are  the  direct  derivatives  Cm^,  Cm^,  Cl^,  Cl^,  Cfp,  Cn^,  and  Cy^;  the  cross  derivatives 
Cnp,  and  Cyp?  and  the  cross-coupling  derivatives  Cf^,  C„^,  Cmp,  and  The  missile 
configurations  used  for  the  investigation  are  representative  of  the  conventional  yaw-to-turn 
and  bank-to-turn  missiles.  The  sensitivity  of  each  configuration  is  documented  with  respect 
to  derivative  variations  Root  loci,  time  response,  and  damping  ratio  plots  are  used. 

2.0  METHOD  OF  ANALYSIS 

The  sensitivity  of  flight  vehicle  motion  to  variations  in  aerodynamic  and  physical 
characteristics  is  reflected  in  the  longitudinal  and  lateral-directional  stability  modes.  The  two 
longitudinal  modes  are  short  period  and  phugoid,  and  the  three  lateral-directional  modes  are 
roll,  spiral,  and  dutch  roll.  Linearized  analysis  of  geometrically  and  aerodynamically 
symmetric  air  vehicles  in  wings-level,  nonrotating  flight  without  sideslip  does  not  require 
more  than  a  three-degree-of-freedom  program  since  the  longitudinal  and  lateral-directional 
perturbations  are  truly  uncoupled.  However,  for  asymmetric  geometry  and  nonzero- 
rotational  rates,  the  perturbations  from  the  reference  conditions  are  all  coupled.  In  many 
cases,  the  degree  of  coupling  may  be  small;  hence,  a  three-degree-of-freedom  analysis  would 
be  adequate,  but  when  an  attempt  is  made  to  determine  the  actual  effect  of  the  cross¬ 
coupling  between  the  longitudinal  and  lateral-directional  motions,  a  six-degree-of-freedom 
analysis  is  required. 

The  Arbitrary  Degree  of  Freedom  (ADOF)  computer  program  was  modified  to  provide 
initial  trim  required  for  level  and  steady  turning  flight  conditions.  The  program  is  applicable 
to  flight  vehicles  either  with  or  without  feedback  control  systems  and  can  be  used  for  any  set 
of  degrees  of  freedom.  The  equations  of  motion  and  data  input  routines  used  in  the  program 
are  of  a  general  nature  and  are  not  restricted  by  assuming  a  symmetric  vehicle  or  an 
equilibrium  reference  condition. 
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The  modes  of  motion  are  represented  by  the  roots  of  the  characteristic  equation  formed 
by  linearization  of  the  six  equations  of  motion. 

Initially,  all  dynamic  derivatives  under  investigation  were  individually  assigned  a  value  of 
unity  and  then  varied  over  a  predetermined  range.  All  other  pertinent  parameters  were 
maintained  at  the  required  trim  values.  The  stability  characteristics  were  recalculated  for 
each  assigned  value  of  the  derivative,  thereby  mapping  the  effects  on  the  five  primary  modes 
of  motion. 


3,0  TECHNICAL  DATA 


3.1  MISSILE  CONFIGURATIONS 

Two  missile  configurations  were  selected  for  the  dynamic  derivative  sensitivity  study  —  a 
conventional  yaw-to-turn  and  a  bank-to-tum  configuration.  The  selection  of  a  specific 
missile  for  each  configuration  was  made  on  the  basis  of  current  configuration  design  trends 
and  available  wind  tunnel  data. 

For  the  bank-to-turn  missile,  one  of  the  early  Interlab  Air-to-Air  Technology  (ILAAT) 
designs  was  selected  and  is  shown  in  Fig.  1,  along  with  representative  full-scale  dimensional 
characteristics.  The  basic  configurational  components  consisted  of  nose,  body,  wing,  and 
fins.  The  body  makes  a  smooth  transition  from  circular  to  an  increasing  elliptic  cross  section 
to  a  constant  elliptical  main  body. 

A  general  research  model  known  as  the  Aerodynamic  Data  Correlation  (ADC)  missile 
was  selected  for  the  yaw-to-turn  configuration.  Representative  full-scale  dimensional 
characteristics  are  given  in  Fig.  2.  The  ADC  model  consists  of  an  ogive  nose  section  and  a 
cylindrical  body  with  aft  mounted  cruciform  fins. 

3.2  AERODYNAMIC  DATA 
3.2.1  General 

The  rigid-body  aerodynamic  data  are  input  to  the  stability  axis  system.  Data  are  given  in 
table  look-up  form  as  functions  of  the  variables  shown  in  the  equations  of  the  Appendix.  All 
aerodynamic  coefficients  for  both  missile  configurations  are  referenced  to  body  cross- 
sectional  area  and  diameter. 
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3.2.2  M^ind  Tunnel  Measurements 

The  static  aerodynamic  data  matrices  used  in  modeling  both  missile  configurations  were 
obtained  from  wind  tunnel  tests  at  AEDC.  The  static  data  for  the  bank-to-turn 
configuration  (Fig.  1)  were  obtained  for  the  Mach  number  range  from  0,8  to  3.5  and  for  an 
angle  of  attack  of  up  to  26  deg.  The  static  data  matrices  for  the  yaw-to-turn  configuration 
(Fig.  2)  were  formulated  from  data  obtained  in  previous  study  at  AEDC,  The  matrices 
include  data  for  Mach  numbers  from  0.6  to  3.0  and  for  an  angle  of  attack  from  -6  through 
26  deg. 

The  matrices  for  the  dynamic  aerodynamic  data  for  the  yaw-to-turn  configuration  were 
formulated  from  data  presented  in  Refs.  1  and  2  and  in  other  work  at  AEDC.  These 
matrices  include  data  in  the  Mach  number  range  from  0.6  to  3.0  and  in  the  angle-of-attack 
range  from  -6  through  26  deg.  The  roll  damping  characteristics  are  documented  in  Refs.  1 
and  2,  and  pitch  and  yaw  damping  have  been  documented. 

3.2.3  Estimated 

The  dynamic  stability  characteristics  for  the  bank-to-turn  configuration  have  not  been 
determined  experimentally.  Several  prediction  methods  were  reviewed  for  estimating  flight 
vehicle  dynamic  derivatives.  Both  semi-empirical  (Ref,  3)  and  unsteady  panel  methods  (Ref. 
4)  were  reviewed.  For  the  methods  available,  a  semi-empirical  method  known  as  the  USAF 
Stability  and  Control  Datcom  (Ref.  3)  appeared  to  be  the  most  extensive  and  easiest  to 
implement.  A  portion  of  Datcom  (Digital  Datcom;  see  Ref.  5)’was  used  in  estimating  the 
bank-to-turn  dynamic  stability  derivatives.  The  Mach  number  and  angle-of-attack  ranges 
for  which  the  derivatives  were  estimated  were,  respectively,  from  0.8  to  3.5  and  from  -6  to  26 
deg.  It  should  be  noted  that  Datcom  is  valid  only  for  the  attached  flow  regime  (low  angle  of 
attack). 

The  missile  configuration  of  Fig.  1  was  geometrically  modified  to  conform  to  the 
restraints  of  the  Digital  Datcom  computer  code.  In  the  Digital  Datcom,  flight  vehicles  are 
modeled  with  the  use  of  body,  wing,  horizontal,  and  vertical  tail  input  data.  Therefore,  for 
the  computer  program,  fins  3  and  4  of  the  missile  of  Fig.  1  were  modeled  as  the  horizontal 
control  surfaces,  while  the  areas  of  fins  1  and  2  (also  of  Fig.  1),  were  combined  as  one  upper 
vertical  surface.  In  addition,  the  wing,  horizontal,  and  vertical  tail  were  modeled  as  straight 
tapered  planforms  rather  than  double  delta  planforms. 
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4.0  RESULTS  AND  DISCUSSION 


4.1  GENERAL 

The  dynamic  stability  analyses  for  the  bank-to-turn  and  yaw-to-turn  missile 
configurations  were  conducted  in  trimmed  level  and  steady  turning  flight.  Turning  flight 
conditions  were  used  to  simulate  the  missiles  in  maneuvering  flight  regime.  The  turning 
flight  condition  provides  larger  values  of  the  angular  rate  terms  (p,  q,  and  r),  thereby  giving 
the  aerodynamic  cross-coupling  moments  larger  values  in  the  linearized  equations  of 
motion.  This  turning  flight  condition  provides  a  means  for  ascertaining  the  importance  of 
the  second-order  cross-coupling  derivative  terms. 

The  dynamic  derivative  sensitivity  studies  were  conducted  with  use  of  flight  envelopes 
representative  of  both  missile  categories.  The  specific  Mach  number,  altitude,  and  g  flight 
conditions  at  which  the  investigations  were  conducted  are  shown  in  Fig.  3.  The  trimmed 
flight  conditions  shown  for  each  Mach  number  and  altitude  are  1-g-level,  an  intermediate  g, 
and  a  maximum  g.  The  maximum  g  at  each  Mach  number  and  altitude  was  limited  only  by 
the  control  surface  authority  and  by  the  maximum  angle  of  attack  of  the  data  matrix.  Tables 
1  and  2  summarize  the  aerodynamics  used  for  the  missile  flight  conditions.  Table  3  gives  the 
mass,  inertia,  and  geometric  characteristics  used  to  represent  the  full-scale  bank-to-turn  and 
yaw-to-turn  missiles. 

Figure  4  gives  the  range  for  varying  the  dynamic  derivatives  plus  the  associated  maximum 
and  minimum  values  measured  in  the  wind  tunnel  and  estimated  by  Digital  Datcom.  Both 
the  direct  and  cross  derivative  ranges  were  arbitrarily  selected  at  approximately  twice  the 
maximum  of  the  measured  and  estimated  derivatives.  The  range  for  the  cross-coupling 
derivatives  was  chosen  to  correspond  to  the  range  for  the  cross  derivatives.  Experimental 
results  obtained  by  Orlik-Ruckemann  (Ref.  6)  with  use  of  a  cone  wing  model  have  shown 
that  the  cross-coupling  derivatives  combined  with  the  acceleration  derivatives  (Cn,^,  Cn^, 
and  Cy)  approach  and/or  exceed  those  of  the  cross  derivatives  for  high  angles  of  attack.  It 
should  be  noted  that  the  dynamic  derivatives  for  the  yaw-to-turn  configuration  were 
experimentally  obtained  as  a  combination  of  the  rate  and  acceleration  terms  (e.g.,  Cn^  -  Cn^ 
cos  a).  In  the  sensitivity  study,  these  combination  derivative  values  were  used  as  pure  rate 
terms. 

The  root  locus  format  shown  in  Fig.  5  demonstrates  the  primary  method  used  to  show 
the  effect  of  derivative  variations  on  the  vehicle  modes  of  motion.  The  resultant  roots 
(eigenvalues)  of  the  characteristic  equation  are  plotted  in  the  root  locus  format  shown  as  the 
derivatives  are  varied.  The  real  part  of  the  root  represents  the  damping  ratio,  whereas  the 
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imaginary  part  represents  the  frequency  of  oscillation  (radian /sec).  The  time  required  for 
each  mode  to  damp  to  one  half  or  to  diverge  to  double  amplitude  is  shown  below  each  plot; 
the  motion  period  is  shown  vertically.  The  roots  of  the  characteristic  equation  may  be  either 
real  (aperiodic)  or  complex  (oscillatory)  and  either  stable  or  unstable. 

4.2  BANK-TO-TURN  CONFIGURATION 

The  dynamic  derivative  sensitivity  analysis  was  conducted  for  altitudes  of  10,000  and 
40,000  ft  for  Mach  numbers  of  0,8.,  1.1,  and  3.5  for  load  factors  of  up  to  25  g.  The  low 
Mach  number  flight  conditions  correspond  to  initial  maneuvers  after  launch,  whereas  the 
high  Mach  number  corresponds  to  supersonic  terminal  flight  conditions. 

4.2.1  Dynamic  Direct  Derivative  Variations 

Longitudinai.  The  missile  motion  sensitivity  to  variations  in  the  direct  rate  and  acceleration 
derivatives  and  is  shown  in  Figs.  6  and  7,  respectively.  The  range  (±  1,000  per 
radian)  for  which  the  derivatives  were  varied  was  at  least  twice  the  magnitude  of  that 
obtained  from  Digital  Datcom  (see  Fig.  4).  The  influences  of  altitude  and  load  factor  (g’s) 
on  the  missile’s  sensitivity  to  derivative  variations  are  presented  for  each  Mach  number.  As 
previously  noted,  the  maximum  g’s  available  at  each  Mach  number  are  limited  by  the  control 
surface  authority  and  angle^of-attack  limit.  Only  the  short  period,  S/P,  mode  shows 
sensitivity  to  variations  in  the  pitching  moment  caused  by  the  and  derivatives.  The 
sensitivity  in  the  short  period  mode  at  each  Mach  number  is  essentially  unaffected  by  the 
load  factor,  g;  there  is,  however,  a  significant  effect  of  Mach  number  on  sensitivity  of  the 
short  period  mode  with  and  variations.  The  short  period  mode  is  less  sensitive  to 
Cn,^  and  Cmi  variations  at  40,000  ft  than  at  10,000  ft.  The  reduced  sensitivity  is  primarily 
due  to  the  reduction  in  density  with  increased  altitude. 

The  motion  sensitivity  to  the  force  derivatives  Ctq  and  is  shown  in  Fig.  8.  The 
longitudinal  and  lateral-directional  modes  of  motion  were  totally  insensitive  to  variations  in 
either  CLq  or  Cl^  over  the  range  of  ±  500  per  radian. 

Figures  9  and  10  demonstrate,  respectively,  the  influence  of  variations  in  Cmq  and 
on  the  vehicle  damping  ratio  and  motion  response  characteristics  at  Mach  3.5  and  10,000-ft 
altitude  for  a  10-g  banked  turn.  Figures  9a  and  10a  show  the  sensitivity  of  the  short  period 
and  damping  ratio,  Z,  with  variations  in  Cmq  and  Cm^  over  the  range  of  ±  1,000  per  radian. 
As  can  be  observed  in  the  root  locus  plots,  the  short  period  damping  ratio  is  quite  sensitive 
to  variations  in  Cmq  and  Cm^.  The  time  response  shown  in  figures  9b  and  10b  results  from  a 
pitch  control  doublet  and  characterizes  the  changes  in  the  short  period  motion  for  three 
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values  each  for  both  Cm^  and  Cfn<j.  The  changes  in  the  pitch  rate  response  to  a  unit  doublet 
in  pitch  control  are  shown  for  values  of  Cm^  and  Cm^  corresponding  to  tlie  trim  values, 
values  near  the  point  at  which  the  short  period  is  neutrally  stable,  and  values  near  the 
maximum  or  minimum  for  the  experimental  ranges  shown  in  Fig.  4.  The  changes  in  the  time 
response,  again,  indicate  the  senstivity  of  the  short  period  motion  to  the  direct  derivatives 
Cn,q  and  Ctn^  and  quantify  the  significance  of  the  root  migration  with  derivative  variation 
shown  in  Figs.  6  and  7.  , 

Lateral-Directional.  The  missile  motion  sensitivity  to  variations  (±1,000  per  radian)  in  the 
direct  damping  derivative  in  yaw,  Cn^,  is  shown  in'Fig.  11.  The  only  significant  motion 
senstivity  to  Cnj.  is  in  the  dutch  roll  mode  of  motion.  The  sensitivity  in  the  dutch  roll  mode 
due  to  increases  with  increasing  load  factor  (g’s)  for  all  Mach  numbers  at  each  altitude. 

The  effects  of  variations  (±  1,0(X)  per  radian)  in  the  direct  roll  damping  derivative,  C(^, 
on  the  longitudinal  and  lateral-directional  modes  of  motion  are  shown  in  Fig.  12,  Only  the 
lateral-directional  modes  are  shown  here  for  clarity;  the  longitudinal  modes  were  insensitive 
to  Cfp  variations. 

The  extreme  sensitivity  of  both  the  roll  and  dutch  roll  modes  to  variations  in  Cf  is 
associated  with  the  low  moment  of  inertia,  I*,  about  the  roll  axis.  The  sensitivity  in  the  dutch 
roll  mode  due  to  Cfp  variations  increases  with  increasing  load  factor  (g’s)  for  all  Mach 
numbers  at  each  altitude.  For  large  negative  values  of  Qp,  the  dutch  roll  mode  may  migrate 
from  a  periodic  spiral  to  a  roll-spiral  coupling  (lateral  phugoid)  to,  finally,  an  aperiodic 
mode  at  large  positive  values.  This  sensitivity  is  made  apparent  by  tracing  the  dutch  roll  root 
migration  at  Mach  0.8  and  40,000-ft  altitude  for  a  5-g  banked  turn.  The  roll  mode  is 
extremely  sensitive  to  relatively  small  changes  in  Cfp  as  can  be  seen  by  tracing  the  roll  root 
migration  at  M  =  3,5  and  altitude  =  10,000  ft  for  1-g  level  flight.  The  sensitivity  in  the  roll 
mode  due  to  Cfp  increases  with  increasing  Mach  number  but  decreases  with  increasing 
altitude  and  load  factor.  The  roll  and  dutch  roll  modes  are  more  sensitive  to  the  same 
magnitude  variations  in  Cfp  than  to  variations  in  Cn^;  this  greater  sensitivity  to  Cfp  can  be 
seen  by  comparing  Figs.  1 1  and  12.  The  increased  sensitivity  of  the  roll  and  dutch  roil  modes 
to  Cfp  is  associated  with  the  previously  mentioned  low  moment  of  inertia,  1*,  about  the  roll 
axis. 

Figure  8  shows  the  motion  sensitivity  to  variations  in  the  side  force  derivative,  Cy^,  with 
variations  in  Ci^  and  Ci^.  The  longitudinal  and  lateral-directional  modes  of  motion  were 
totally  insensitive  to  Cy^  over  the  range  of  ±  500  per  radian. 

The  influence  of  variations  in  Cn^  on  the  vehicle  lateral-directional  damping  ratio  and 
motion  response  is  shown  in  Fig.  13  at  M  =  0.8  and  40,000-ft  altitude  for  a  10-g  banked 
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turn.  The  sensitivity  of  the  dutch  roll  and  spiral  mode  damping  ratio  to  variations  in  C„^  is 
'  shown  in  Fig.  13a.  As  can  be  noted  in  the  root  locus  plots  of  Fig.  1 1 ,  the  dutch  roll  damping 
is  very  sensitive  to  Cn,  variations.  The  yaw  rate  time  response  shown  in  Fig.  13b  results  from 
a  yaw  control  doublet  and  characterizes  the  changes  in  the  dutch  roll  motion  for  three  values 
of  Cn^.  The  sensitivity  of  the  dutch  roll  mode,  as  demonstrated  by  the  yaw  rate  response, 
ranges  from  a  highly  damped  motion  for  Cn^  =  -1,000  per  radian  to  a  lightly  damped  motion 
that  approaches  a  neutrally  stable  mode  at  Cn^  =  100  per  radian.  The  change  in  the  yaw  rate 
response  for  the  three  values  of  C„,  corresponds  to  the  changes  in  the  Time-to-Half- 
Amplitude,  T|/2,  shown  in  Fig.  11  for  the  same  values  of  Cn^. 

The  sensitivity  of  the  dutch  roll  mode  damping  ratio  and  motion  response  to  variations  in 
Qj,  over  the  range  of  ±1,000  per  radian  is  presented  in  Fig.  14.  The  flight  condition 
presented  corresponds  to  M  =  0.8  and  altitude  =  10,000  ft  for  a  10-g  banked  turn.  The 
damping  ratio  plot  (Fig.  14a)  demonstrates  the  sensitivity  of  the  dutch  roll  to  relatively  small 
changes  in  Cfp  about  zero  value.  This  sensitivity  is  also  apparent  in  the  roll  rate  time  response 
(Fig.  14b)  for  values  of  -200,  -100,  and  -62.5  (trim  value)  per  radian.  The  roll  rate 
response  attributable  to  a  unit  roll  control  doublet  changes  from  a  motion  with  T1/2  ~  0.1 
sec  at  Cfp  =  -62.5  to  a  much  less  damped  motion  with  T1/2  -  0.04  sec  at  =  -200  per 
radian.  The  roll  and  dutch  roll  modes  are  quite  sensitive  to  small  changes  in  Qp^  this  can  be 
seen  by  tracing  the  root  migrations  in  Fig.  12a. 

Both  the  roll  and  dutch  roll  modes  of  motion  have  been  shown  to  be  quite  sensitive  to 
changes  in  and  Cfp  within  the  value  ranges  estimated  by  Datcom  (Fig.  4). 

4.2.2  Dynamic  Cross  Derivative  Variations 

The  missile  motion  sensitivity  to  variations  in  the  cross  derivative,  Ct^,  is  shown  in  Fig. 
15.  Only  the  lateral-directional  modes  of  motion  are  affected  by  a  variation  in  (in  the 
range  of  ±500  per  radian).  The  sensitivity  of  the  dutch  roll  mode  to  increases  with 
increasing  load  factor  (g’s)  with  each  Mach  number,  whereas  the  roll  mode  sensitivity 
decreases  with  increasing  g  flight  at  each  altitude  and  Mach  number.  By  comparing  the 
dutch  roll  root  migrations  in  Figs.  11  and  15,  one  sees  that  the  sensitivity  of  the  dutch  roll 
mode  to  variations  in  Cf^  is  equal  to  that  for  variations  in  the  direct  derivative,  Cnp- 

Variations  of  the  yawing  moment  caused  by  the  roll  rate  parameter,  C„p,  affect  the 
damping  and  frequency  of  the  roll  and  dutch  roll  modes  of  motion  significantly,  as  can  be 
seen  by  tracing  the  root  migration  in  Fig,  16.  In  general,  the  sensitivity  of  the  roll  and  dutch 
roil  modes  to  C^p  increases  with  increasing  load  factor  (g’s)  and  Mach  number.  The  roll 
mode  may  go  from  aperiodic  stable  with  negative  values  of  Cnp  to  unstable  with  positive 
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values  of  C„p,  The  dutch  roll  mode  may  go  from  an  unstable  oscillation  with  negative  values 
of  Cnp  to  a  stable  oscillation  and  finally  to  a  degenerated  aperiodic  mode  with  positive  values 
of  Cnp.  In  general,  the  roll  and  dutch  roll  modes  of  motion  are  more  sensitive  to  C,ip  (Fig.  16) 
than  to  C(y  (Fig.  15)  for  the  flight  conditions  investigated. 

The  motion  sensitivity  to  variations  in  the  side  force  derivative,  Cyp,  is  shown  in  Fig.  17. 
As  with  other  force  derivatives,  the  effect  of  Cyp  on  both  the  roll  and  dutch  roll  modes  is  not 
considered  significant. 

The  influence  of  variations  in  Ci,.  on  the  vehicle  dutch  roll  damping  ratio  and  on  the  roll 
rate  response  is  shown  in  Fig.  18  at  M  =  0.8  and  altitude  =  40,000  ft  for  a  10-g-banked 
turn.  Shown  in  Fig.  18a  is  the  sensitivity  of  the  dutch  roll  mode  damping  ratio  to  variations 
in  Cf^  over  the  range  ±  500  per  radian.  The  roll  rate  time  responses  shown  in  Fig.  18b  result 
from  a  yaw  control  doublet  and  demonstrate  the  sensitivity  of  the  dutch  roll  damping  to 
C;^.  The  dutch  roll  motion,  as  shown  by  the  roll  rate  response,  ranges  from  highly  damped 
motion  for  C(y  =  500  to  that  approaching  a  neutrally  stable  mode  at  Ct,.  =  -50  per  radian. 
The  changes  in  the  roll  rate  time  responses,  along  with  the  dutch  roll  root  migrations  shown 
in  Fig.  15,  indicate  that  the  missile’s  motion  may  be  highly  sensitive  to  variations  in  the  cross 
derivative  C,'^. 

The  sensitivity  of  the  dutch  roll  damping  ratio  and  yaw  rate  response  to  Cnp  is  presented 
in  Fig.  19  at  M  =  1.1  and  altitude  =  10,000  ft  for  a  10-g  banked  turn.  The  dutch  roll  damp¬ 
ing  ratio  (Fig.  19a)  is  sensitive  to  variations  in  Cnp  over  the  range  of  ±  500  per  radian.  The 
yaw  rate  time  response  shown  in  Fig.  19b  results  from  a  roll  control  double!  and 
demonstrates  the  sensitivity  of  the  dutch  roll  mode  motion  for  three  values  of  Cnp.  The 
changes  in  the  yaw  rate  time  response  with  Cpp  variations  follows  trends  of  the  dutch  roll 
mode  root  migrations  with  Cnp  variations  as  shown  in  Fig.  16.  As  with  Qy  the  bank-to-turn 
missile  motion  is  sensitive  to  Cnp  variations. 

The  influence  of  variations  in  Cyp  on  the  dutch  roll  damping  ratio  and  yaw  rate  response 
is  shown  in  Fig.  20  at  M  =  3.5  and  altitude  =  10,000  ft  for  a  25-g  banked  turn.  For  variation 
in  Cyp  over  the  range  of  ±  500  per  radian,  the  change  in  the  dutch  roll  damping  ratio  of  Fig. 
20a  is  not  considered  large  enough  to  significantly  alter  the  dutch  roll  motion.  This  is 
shown  in  the  minimal  change  in  the  yaw  rate  response  of  Fig.  20b  attributable  to  a  roll 
control  doublet  for  range  in  Cyp  of  -500  to  500  per  radian. 

In  general,  therefore,  direct  damping  derivatives  are  important  in  motion  analysis.  The 
sensitivity  of  the  missile  lateral-directional  modes  of  motion  to  variations  in  the  cross 
derivatives  Ci,.  and  Cnp  is  on  the  order  of  the  direct  derivative  Cn,;  therefore,  and  C^p  are 
necessary  to  perform  accurate  missile  motion  simulation. 
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4.2.3  Dynamic  Cross>Coupling  Derivative  Variations 

In  the  cross-coupling  derivative  investigation  the  derivatives  Cm^,  Cmp,  ,  and  were 
each  initially  given  nominal  values  of  100  per  radian.  Then,  individual  variations  in  each 
derivative  were  performed  with  the  remaining  derivatives  at  their  nominal  values.  By 
maintaining  nominal  (nonzero)  values  for  the  cross-coupling  derivatives,  the  possible 
derivative  interaction  effects  described  in  Ref.  7  were  included. 

The  motion  sensitivity  to  variations  in  the  cross-coupling  derivative  is  shown  in  Fig. 
21  for  Mach  numbers  of  0.8,  1.1,  and  3.5.  Including  the  cross-coupling  derivatives  in  the 
linearized  analysis  results  in  a  weak  coupling  in  the  longitudinal  and  lateral-directional  axes 
as  indicated  by  the  migration  in  the  short  period,  spiral,  roll,  and  dutch  roll  modes.  As 
observed,  the  sensitivity  of  the  dutch  roll  mode  with  variations  in  depends  on  the  flight 
conditions  (Mach  number,  altitude,  and  load  factor).  Although  the  sensitivity  of  the  modes 
of  motion  to  is  not  as  significant  as  it  is  to  (compare  Figs.  15  and  21),  errors  may 
occur  in  missile  motion  simulation  if  the  derivative  is  not  modeled  when  the  magnitude 
approaches  the  maximum  ranges  utilized  in  Fig.  21. 

The  motion  sensitivity  with  variations  in  the  cross-coupling  derivative  Cmp  is  shown  in 
Fig.  22.  The  short  period,  roll,  and  dutch  roll  mode  sensitivity  to  C^p  is  similar  to  that 
experienced  with  the  derivative  variations.  If  it  is  assumed  that  the  cross-coupling 
derivative  C^p  is  of  a  magnitude  approaching  that  of  the  cross  derivatives,  then  incorrect 
motion  simulation  may  possibly  result  if  the  derivative  is  not  included  in  the  aerodynamic 
modeling. 

The  motion  sensitivity  with  variations  in  the  cross-coupling  derivatives  and  is 
shown  in  Fig,  8  along  with  the  variations  in  the  force  derivatives.  The  longitudinal  and 
lateral-directional  modes  of  motion  are  totally  insensitive  to  variations  in  either  Cm^  or  Cp 
over  the  range  of  ±  500  per  radian. 

The  sensitivities  of  the  missile  damping  ratio  and  frequency  and  the  roll  and  pitch  rate  to 
variations  of  (Fig.  23)  and  Cmp  (Fig,  24)  are  given  for  M  =  3.5  and  altitude  =  40,000  ft 
for  a  25-g  banked  turn.  At  this  flight  condition,  the  missile  was  most  sensitive  to  variations 
in  cross-coupling  derivatives.  The  short  period,  phugoid,  and  dutch  roll  mode  damping 
ratios  are  shown,  in  Fig.  23a,  to  be  quite  sensitive  to  variations  in  Cf^^  over  the  range  of  ±  500 
per  radian.  The  short  period  and  dutch  roll  damping  ratios  exhibit  similar  sensitivity  to  Cmp 
(Fig.  24b).  The  sensitivity  of  the  dutch  roll  mode  of  motion  to  C^  is  shown  (Fig.  23b)  by  the 
changes  in  the  roll  rate  response  resulting  from  a  pitch  control  doublet,  whereas  the 
sensitivity  of  the  short  period  mode  of  motion  to  Cmp  is  shown  (Fig.  25b)  by  the  changes  in 
the  pitch  rate  response  resulting  from  a  roll  control  doublet. 
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4.3  YAW-TO-TURN  CONFIGURATION 

The  dynamic  sensitivity  study  for  the  yaw-to-turn  missile  was  conducted  at  altitudes  of 
10,000  and  20,000  ft  at  Mach  1 .3  and  at  altitudes  of  10,000  and  40,000  ft  for  Mach  3.0.  The 
flight  conditions  investigated  were  1-g-level  flight  at  Mach  1.3  and  1-g-Ievel  and  5-  and  10-g 
turning  flight  at  Mach  3.0.  Angle  of  attack  and/or  control  surface  deflection  required  for 
trim  flight  at  Mach  numbers  of  less  than  1 .3  for  accelerated  g  flight  conditions  exceeded  the 
range  in  the  available  aerodynamic  data  matrix. 

4.3.1  Dynamic  Direct  Derivative  Variations 

Longitudinal.  Figures  25  and  26  show  the  missile  motion  sensitivity  to  variations  in  the  direct 
rate  and  acceleration  derivatives  Cmg  and  Cm^.  The  range  for  which  the  derivatives  were 
varied  (+1,000  per  radian)  was  slightly  less  than  twice  the  magnitude  obtained 
experimentally  (see  Fig.  4).  Due  to  the  absence  of  coupling  between  the  longitudinal  and 
lateral-directional  planes  of  motion,  only  the  longitudinal  modes  are  shown.  Only  the  short 
period  mode  shows  sensitivity  to  variations  in  Cm^  and  Cm^.  There  is  a  significant  Mach 
number  effect  on  the  sensitivity  of  the  short  period  mode  with  and  Cm^  variations.  As 
noted  for  the  bank-to-turn  configuration,  the  short  period  mode  is  less  sensitive  at  40,000  ft 
to  Cmq  and  Cm^  variations  than  at  10,000  ft  because  of  the  decrease  in  density  with 
increasing  altitude. 

The  motion  sensitivity  to  the  force  derivatives  Ci^  and  Cl^  is  shown  in  Fig.  27.  The 
missile  modes  of  motion  were  totally  insensitive  to  variations  in  either  Cl^  or  Cl^  over  the 
range  of  ±  500  per  radian. 

The  influence  of  variations  in  Cm,  and  Cm^  on  the  vehicle  longitudinal  damping  ratio 
and  pitch  rate  response  is  shown  respectively  in  Figs.  28  and  29  at  M  =  3.0  and  altitude  = 
40,000  ft  for  1-g-level  flight.  In  Figs.  28a  and  29a,  the  short  period  damping  ratio  is  shown  to 
be  sensitive  to  variations  in  and  The  pitch  rate  time  response  shown  in  these 
figures  results  from  a  pitch  control  doublet  and  demonstrates  the  sensitivity  of  the  short 
period  mode  to  and  Cm^.  The  changes  shown  in  the  pitch  rate  response  correspond  to 
values  of  and  Cm^  at  trim,  at  values  of  zero,  and  at  values  near  which  the  short  period 
mode  is  neutrally  stable. 

Lateral-Directional.  The  missile  motion  sensitivity  to  variations  (±  1,0(K}  per  radian)  in  the 
direct  damping  derivatives  in  yaw,  C„^,  is  shown  in  Fig.  30.  The  roll,  spiral,  and  dutch  roll 
modes  are  affected  to  different  degrees  by  Cn^  variation.  The  most  predominant  effect  is  on 
the  dutch  roll  mode  for  the  flight  conditions  evaluated.  The  sensitivity  of  the  roll  and  spiral 
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modes  is  more  pronounced  for  the  1-g-level  trim  flight  conditions.  The  dutch  roll  mode  sen¬ 
sitivity  due  to  Cn^  variations  increases  with  increasing  load  factor  (g’s)  and  altitude  at  Mach 
number  3.0. 

The  effect  of  variations  ( ±  1,000  per  radian)  in  the  direct  roll  damping  derivative,  on 
the  lateral-directional  modes  of  motion  is  shown  in  Fig.  31.  The  extreme  sensitivity  of  the 
roll  and  spiral  modes  to  is  associated  with  the  low  moment  of  inertia,  lx,  about  the  roll 
axis.  The  pitch,  ly,  and  yaw,  Iz,  moment  of  inertia  are  approximately  two  orders  of 
magnitude  larger  (see  Table  3)  than  the  rolling  moment  of  inertia.  The  dutch  roll  mode  of 
motion  is  sensitive  to  variations  in  but  this  sensitivity  does  not  produce  an  unstable 
mode  for  the  flight  conditions  analyzed.  The  roll  and  spiral  modes  are  more  sensitive  to  the 
same  magnitude  variations  in  than  to  Cn^  as  can  be  seen  by  tracing  the  root  migrations  in 
Figs.  30  and  31. 

The  lateral-directional  modes  of  motion  are  totally  insensitive  to  variations  (±500  per 
radian)  in  the  side  force  derivative  Cy^;  these  modes  are  shown  with  the  normal  force 
derivatives  Ci^  and  Cl,;  in  Fig.  27. 

The  influence  of  variations  in  C„^  on  the  vehicle  spiral  and  dutch  roll  damping  and  yaw 
rate  response  is  shown  in  Fig.  32  at  M  =  3.0  and  altitude  =  10,000  ft  for  a  10-g  banked  turn. 
The  sensitivity  of  the  dutch  roll  damping  ratio  with  variations  ( ±  1 ,000  per  radian)  is  shown 
in  Fig.  32a.  As  can  be  noted  from  the  curve  slope,  the  dutch  roll  damping  is  very  sensitive  to 
Cn,  variations.  The  time  history  shown  in  Fig.  32b  results  from  a  yaw  control  doublet  and 
characterizes  the  effect  of  changes  in  Cn,.  on  the  dutch  roll  motion. 

The  sensitivity  of  the  roll  and  dutch  roll  damping  and  roll  rate  response  to  variations  in 

is  shown  in  Fig.  33  at  M  =  3.0  and  altitude  =  10,000  ft  for  a  10-g  banked  turn.  The 
damping  ratio  plot  (Fig.  33a)  demonstrates  the  sensitivity  of  the  spiral,  roll,  and  dutch  roll 
modes  to  relatively  small  changes  in  Ct^  about  zero  values.  This  sensitivity  in  the  dutch  roll 
mode  is  apparent  from  the  roll  rate  time  response  plot  (Fig.  33b).  For  the  same  magnitude 
change  in  the  derivatives,  the  dutch  roll  mode  of  motion  is  significantly  more  sensitive  to 
than  to  Cn^. 

4.3.2  Dynamic  Cross  Derivative  Variations 

The  vehicle  motion  sensitivity  to  the  cross  derivative  is  shown  in  Fig.  34.  Only  the 
lateral-directional  modes  of  motion  are  affected  by  variations  in  over  the  range  of  ±  500 
per  radian.  The  sensitivity  of  the  dutch  roll  mode  to  C(^  increases  with  increasing  load  factor 
(g’s),  whereas  the  spiral  mode  sensitivity  decreases  with  increasing  load  factor.  By 
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comparing  the  dutch  roll  root  migrations  with  Fig.  30,  it  can  be  seen  that  the  sensitivity  of 
the  dutch  roll  mode  to  is  greater  than  that  for  the  direct  derivative 

Changes  in  the  spiral,  roll,  and  dutch  roll  modes  of  motion  with  variation  of  Cnp  are 
shown  in  Fig.  35.  The  roll  and  spiral  mode  sensitivity  to  Cnp  increases  with  increasing  load 
factor.  The  dutch  roll  mode  sensitivity  is  also  highly  dependent  on  the  load  factor.  In  most 
cases,  the  roll  and  spiral  modes  are  more  sensitive  to  Cnp  than  to  C|>^  as  can  be  seen  by  tracing 
the  root  migrations  in  Figs.  34  and  35. 

The  lateral-directional  modes  of  motion  are  insensitive  to  variations  ( ±  500  per  radian)  in 
the  side  force  parameter  Cyp  as  shown  in  Fig.  27.  For  the  flight  conditions  investigated,  the 
Cyp  parameter  does  not  appear  to  be  an  important  parameter  for  motion  analysis. 

The  influence  of  variations  in  on  the  lateral-directional  damping  ratios  and  roll  rate 
response  is  shown  in  Fig.  36  at  M  =  1.3  and  altitude  =  20,000  ft  for  1-g-leveI  flight.  The 
damping  ratio  plot  demonstrates  the  sensitivity  of  the  spiral,  roll,  and  dutch  roll  modes  to 
variation  in  Cf^  over  the  range  of  ±  500  per  radian.  The  roll  and  spiral  modes  are  quite 
sensitive  to  Cf^  over  the  range  of  ±500  per  radian.  This  is  shown  by  changes  in  the  roots 
from  a  coupled  root  at  large  negative  values  to  real  non-oscillatory  roots  or  modes  for  values 
less  negative  than  -100  per  radian.  The  roll  rate  time  response  attributable  to  a  yaw  control 
doublet  (Fig.  36b)  demonstrates  the  predominance  of  the  spiral  and  dutch  roll  roots  in  the 
Cfj.  range  from  -100  to  200  per  radian. 

The  sensitivity  of  the  lateral-directional  damping  ratios  and  yaw  rate  time  response  to 
Cnp  is  shown  in  Fig.  37  for  M  =  3.0  and  altitude  =  10,000  ft  for  a  10-g  banked  turn. 
Changes  in  the  damping  ratio  of  the  dutch  roll  mode  shown  in  Fig.  37a  further  demonstrate 
the  sensitivity  of  the  dutch  roll  roots  to  relatively  small  changes  in  Cnp  near  zero.  The 
changes  in  the  yaw  rate  time  response  (Fig.  37b)  with  changes  in  Cnp  point  up  the 
signiFtcance  of  the  root  migrations  shown  in  Fig.  35.  The  changes  in  the  time  response  shown 
for  values  of  Cnp  =  -1.2  and  50  demonstrate  the  dutch  roll  sensitivity  to  Cnp,  whereas  the 
changes  shown  for  values  of  -1.2  and  -30  demonstrate  the  spiral  mode  sensitivity. 

The  sensitivity  of  the  missile  lateral-directional  modes  of  motion  to  variations  in  the 
cross  derivatives  Ci^  and  Cnp  is  of  the  same  magnitude  as  for  the  direct  derivatives  and 
therefore  is  important  and  necessary  for  accurate  motion  simulation. 

4.3.3  Dynamic  Cross-Coupling  Derivative  Variations 

As  noted  for  the  bank-to-turn  configuration,  the  cross-coupling  derivatives  Cn,^,  Cn,p, 
C(^  ,  and  Cn^  were  assigned  nominal  values  of  100  per  radian.  Each  derivative  was  varied 
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separately  over  the  range  of  ±  500  per  radian  while  the  others  were  fixed  at  their  nominal 
values. 

The  sensitivity  of  the  longitudinal  and  lateral-directional  modes  of  motion  to  variations 
in  the  roll  due  to  pitch  rate  derivative,  C^,  is  shown  in  Fig.  38.  The  short  period,  dutch  roll, 
roll,  and  spiral  modes  are  shown  to  be  quite  sensitive  to  over  the  range  ±  500  per  radian 
for  M  =  3.0  flight  conditions.  The  dutch  roll  mode  sensitivity  to  increases  with 
increasing  load  factor,  whereas  the  roll  sensitivity  remains  approximately  the  same  with 
increasing  load  factor.  Although  the  lateral-directional  modes  of  motion  are  less  sensitive  to 

than  to  Cf^,  inaccurate  motion  simulation  may  occur  if  the  C(^  derivative  is  of  a  large 
magnitude  and  is  not  included. 

The  motion  sensitivity  with  variations  in  the  pitch  attributable  to  roll  rate  derivative, 
Cmp,  is  shown  in  Fig.  39.  The  mode  sensitivity  to  Cmp  is  almost  identical  to  that  shown  for 
C(^.  Neglecting  C^p  in  motion  analysis  could  cause  inaccuracies  if  the  derivative  is  of  the 
magnitudes  over  the  ranges  investigated  in  this  study. 

The  sensitivity  of  the  missile  modes  of  motion  to  variation  in  Cn^  is  shown  in  Fig.  40.  The 
root  locus  plot  for  M  =  1.3  flight  conditions  is  omitted  because  of  the  insensitivity  of  any  of 
the  modes  to  Cn^.  For  the  M  =  3.0  flight  conditions,  one  can  note  some  movement  of  the 
roots  with  Cnq  variation,  but  the  movement  is  much  less  than  that  associated  with  and 
Cmp  variations. 

The  missile  modes  of  motion  are  totally  insensitive  to  variations  in  and  are  included 
in  Fig.  27  with  the  force  derivatives. 

The  sensitivity  of  the  missile  damping  ratio  and  roll  rate  time  response  to  is  shown  in 
Fig.  41  at  M  =  3.0  and  attitude  =  10,000  ft  for  a  10-g  banked  turn.  This  flight  condition 
represents  the  condition  for  which  the  mode  sensitivity  to  was  the  most  pronounced.  The 
short  period,  phugoid,  dutch  roll,  and  spiral  damping  ratios  in  Fig.  41a  are  shown  to  be 
quite  sensitive  to  variations  in  over  the  range  of  ±  500  per  radian.  The  roll  rate  time 
response  (Fig.  41b)  to  a  pitch  control  doublet  demonstrates  the  sensitivity  of  the  dutch  roll 
mode  to  variations  in  over  the  range  from  -50  to  100  per  radian. 

The  sensitivity  of  the  missile  damping  ratio  and  pitch  rate  time  response  to  Cmp 
variations  over  the  range  ±  500  per  radian  is  shown  in  Fig.  42  at  M  =  3.0  and  altitude  = 
10,000  ft  for  1-g-level  flight.  The  primary  effect  shown  in  Fig.  42a  is  the  sensitivity  of  the 
short  period  damping  ratio.  Associated  with  the  change  in  damping  ratio  is  a  corresponding 
change  in  frequency  as  can  be  seen  by  tracing  the  root  migration  in  Fig.  39.  The  pitch  rate 
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time  response  (Fig.  42b)  demonstrates  the  short  period  sensitivity  to  C^p  values  of  -500,  0, 
and  200  per  radian.  A  point  of  particular  interest  is  the  lack  of  coupling  between  the 
longitudinal  and  lateral-directional  modes  for  Cmp  =  0.0. 

The  sensitivity  of  the  damping  ratio  and  yaw  rate  time  response  to  Cn^  is  shown  in  Fig.  43 
at  M  =  3.0  and  altitude  =  10,000  ft  for  10-g  banked  turn.  Little  change  in  dutch  roll  and 
short  period  damping  ratio  is  seen  in  Fig.  43a.  The  changes  in  the  yaw  rate  time  response 
(Fig.  43b)  for  values  of  Cn^  =  -500,  0,  and  200  demonstrate  the  dutch  roll  sensitivity  to  Cn^. 
The  missile  modes  of  motion  are  much  less  sensitive  to  C„^  than  they  are  to  or  Cmp- 

5.0  CONCLUDING  REMARKS 

Based  on  this  analysis,  the  following  observations  and  conclusions  may  be  made: 

la.  Variations  in  the  longitudinal  direct  derivatives  Cm^  and  significantly 
affect  the  short  period  mode  of  motion  for  both  the  bank-to-turn  and  yaw-to- 
turn  missiles.  The  sensitivity  of  the  short  period  mode  of  motion  with  variations 
in  Cnin  and  Cm-  increases  with  Mach  number  but  is  relatively  unchanged  by 
load  factor  (g’s). 

b.  Variations  in  the  lateral-directional  direct  derivatives  Cn^  and  Cfp  significantly 
alter  the  roll,  spiral,  and  dutch  roll  modes  of  motion  for  both  the  bank-to-turn 
and  yaw-to-turn  missiles.  The  sensitivity  of  the  dutch  roll  mode  increases  with 
increasing  load  factor  (g’s)  and  Mach  number.  The  spiral  and  roll  mode 
sensitivity  may  increase  or  decrease  with  increasing  load  factor  (g’s),  depending 
on  the  missile  configuration. 

2.  Variations  in  the  lateral-directional  cross  derivatives  and  alter  the  bank- 
lo-turn  and  yaw-to-turn  roll,  spiral,  and  dutch  roll  modes  of  motion  in  a 
manner  similar  to  that  resulting  from  variations  of  the  same  magnitude  in  the 
direct  derivatives.  The  roll  and  dutch  roll  motion  sensitivity  increases  with 
increasing  load  factor  (g’s)  and  Mach  number. 

3.  Large  variations  (±500  per  radian)  in  the  cross-coupling  derivatives  C^,  Cn^^, 

^d  Cnip  alter  the  longitudinal  and  lateral-directional  modes  of  motion  for  both 
bank-to-turn  and  yaw-to-tum  missiles.  The  yawing  moment  caused  by  pitch 
rate,  Cn^^,  was  found  to  be  less  important  than  either  or  Cmp.  The  pitching 
moment  associated  with  yaw  derivative,  Cm^,  appears  to  be  insignificant  in 
missile  motion. 
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4.  Varying  the  force  derivatives  Cl^,  Cl^,  Cyp,  and  Cy^  did  not  produce  any 
noticeable  changes  in  the  modes  of  motion  for  either  the  bank-to-turn  or  yaw- 
to-turn  missile. 
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Figure  1.  Bank-to-tum  missile  configuration. 


I 

14 


Fin  3 


DIMENSIONS  IN  FEET 


Figure  2.  Yaw-to-turn  missile  configuration. 
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Figure  3.  Missile  flight  envelope. 


as 


MAXIMUM  OR  MINIMUM 

TABLE  VALUES  CONFIGURATION 


RANGES  OF  VARIABLES 


MEASURED 


YAW-TO-TURN 


-  1000 


-  500 


PER 

RADIAN 


0 


500 


RANGE  ±  1000  ESTIMATED  1  I  BANK-TO-TURN 


Cniqi  Cnr  C^p  Cl^*  0^^  Cnp  Cmp 

Cyp 


Figure  4.  Range  of  derivative  values. 


AEDC-TR.aO-11 


ro 


S  =  <J+  jCJ  (ROOTS  OF  CHARACTERISTIC  EQ.) 
^  IMAG  ( jai) 


RELATIONSHIP  OF  TERMS 
Up 

Z  »  DAMPING  RATIO  =  -  tr/utn  =  o-Zo-^p 
T,^2  «  -0.69315/0- 
P  =  PERIOD  =  2  IT /CD 
jPn(l/Z) 

C  t/2  = 


TIME  TO  HALF,  sec 


TIME  TO  DOUBLE,  sec 

Figure  5.  Sample  root  locus  format. 


AEDC-TR-SO-11 


IMAGINARY  PART  IMAGINARY  PART  IMAGINARY  PART 


AEDC-TR-80-11 


TIME  TO  HALF,  sec  TIME  TO  DOUBLE,  sec 


8.  M  =  0.8 

Figure  6.  Bank-to-turn  configuration  -  locus  of  roots  with  variation. 
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b.  M  =  1.1 
Figure  6.  Continued. 
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Figure  7.  Bank-to-turri  configuration  -  locus  of  roots  with  Cn,,i  variation. 
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Figure  7.  Continued. 
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Figure  10.  Banic-to-turn  configuration  -  effect  of  0^^  on  clamping 
ratio  and  time  response. 
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Figure  11.  Continued. 
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Figure  12.  Continued. 
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Figure  13.  Bank-to-turn  configuration  -  effect  of  C„^  on  damping 
ratio  and  time  response. 
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Figure  14.  Bank-to-turn  configuration  •  effect  of  on  damping 
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with  Cd,  variation. 
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Figure  15.  Continued. 
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Figure  16.  Bank-to-turn  configuration  -  locus  of  roots 
with  C„p  variation. 
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Figure  16.  Continued. 
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Figure  17.  Bank-to-turn  configuration  -  locus  of  roots 
with  Cv_  variation. 
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with  Cmj  variation. 
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Figure  30.  Yaw-to-turn  configuration  -  locus  of  roots 
with  Cnr  variation. 
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Figure  30.  Concluded. 
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Figure  31.  Concluded. 
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Figure  32.  Yaw-to-turn  configuration  •  effect  of  ^  on  damping 
ratio  and  time  response. 
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Figure  33,  Yaw-to-turn  configuration  -  effect  of  on  damping 
ratio. and  time  response. 
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Figure  34.  Yaw-to-turn  configuration  -  locus  of  roots 
with  C£,  variation. 
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Figure  34.  Concluded. 
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Figure  35.  Yaw-to-turn  configuration  -  locus  of  roots 
with  Cnp  variation. 
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Figure  35.  Concluded. 
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Figure  36.  Yaw*to-turn  configuration  -  effect  of  Ci;  on 
damping  ratio  and  time  response. 


85 


AEDC-TR-80-11 


O 


£ 

s 

4 

O 

I 


SYM  WCNITUOt  OF  CAIN 

LARGE  K  GAIN  -  0 

LARGE  &  GAIN  -  1 

LARGE  »  GAIN 


C„p,  pir/rad 

a.  Damping  ratio 


TIME,  tee 

b.  Time  response 

Figure  37.  Yaw-to-turn  configuration  -  effect  of  Cn  on  damping 
ratio  and  time  response. 


86 


IMAGINARY  PART 


12 


8 


0 

-8 


SYM 
lARGE  X 
LARGE  & 
URGE  0 

i 

o 

* 

o 


0.1 


/WAGNITUPE  OF  GAIN 
GAIN  ‘  0 
GAIN  •  1 
GAIN  -  ±  0) 

0<  iLOGifldGAINUt  <1 
1  <  iLOGi^lGAINllt  <2 
2<  iLOGiodGAINllI  <3 

3  <  iLOGiodGAINlll  <4 

4  <  iLOGjiyllGAINlll  <  3 

5  <  ILOGiodGAINDI  <  6 
6<  iLOGiodGAINlH  <  7 
7  <  iLOG^ollGAINDI  <  8 
S  <  ILOGiqIIGAINDI  <  9 


LEVEL  FLIGHT,  Ig 
Qt  i,=6.9d8g 

aI*'-T  =  <  =  9.3d., 

TalT=  20'^ 


D/R, 


\500 

0  ,  S/P, 10 
500 
500 

-500  ^Oi_S/P.  20*' 

500 

500 
-  500 

¥o  D/R,  20*' 

-  [-■ 

PHUGOlO,  lO'',  20^ 


R,  10 
500 

-t 


500 

20 


-^.6 


SPIRAL  ,  10*',  20'‘ 


.75  o 

o 


a 

o 


tr 

UJ 

a. 


5 

10 


-6 


-  4  -2 

REAL  PART 

_ I _ L 


1 


J 


0.25 

TIME  TO  HALF  ,  S«c 


I 


-I  -0.25 

TIME  TO  DOUBLE, 


sec 


a.  M  =  1.3 

Figure  38.  Yaw-to-turn  configuration  -  locus  of  roots 
with  variation. 
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Figure  39.  Yaw-to-tum  configuration  •  locus  of  roots 
with  Cm  p  variation. 
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Figure  43.  Yaw-to-turn  configuration  -  effect  of  Cn,,  on  damping 
ratio  and  time  response. 
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AEDC-TR-80-11 


Table  1.  Bank-To-Turn  Configuration 
a.  M  =  0.8,  Altitude  =  10,000  ft 


Trim  Values 

Level 

Flight,  1  q 

Turning 
Plight,  5  g 

Turning 
Flight,  10 

a (deg) 

1.5 

5.8 

9.8 

iKdeg) 

0 

78.5 

84.3 

deg/sec) 

0 

10.5 

21.3 

{per  rad) 

-264.9 

-286.6 

-301.3 

C  (per  rad) 

0 

0 

0 

Cl 

(per  rad) 

P 

0 

0 

0 

C  (per  rad) 

0 

0 

0 

r 

C,,  (per  rad) 

137.4 

169.8 

186.4 

^ (per  rad) 

0 

0 

0 

a 

C  (per  rad) 

-66.3 

-67.0 

-66.8 

r 

C  (per  rad) 

-1.0 

-7.4 

-13,3 

P 

C  (per  rad) 

0 

0 

0 

q 

C  (per  rad) 

13.4 

13.8 

14.3 

r 

C  (per  rad) 

-62.6 

-62.8 

-62.  5 

C  (per  rad) 

0 

0 

0 

q 

C  (per  rad) 

"r 

0.3 

0.2 

0.2 

C  (per  rad) 

-0.7 

0.6 

1.8 

p 
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Table  1.  Continued 
b.  M  =  0^,  Altitude  =  40,000  ft 


Level  Turning  Turning 


Trim  Values 

Flight,  1  g 

Flight,  5  g 

Flight,  10 

a  <deg} 

4.4 

15.1 

24.7 

ifi  (deg) 

0 

78.5 

84.3 

ji  (deg/sec) 

0 

11.8 

23.9 

C  (per  rad) 

-274.8 

-311.2 

-330.5 

(per  rad) 

0. 

0 

0 

0 

(per  rad) 

0 

0 

0 

C  (per  rad) 

"'r 

0 

0 

0 

Ct  (per  rad) 

157.6 

205.4 

231.0 

Ct  (per  rad) 

0 

0 

0 

C  (per  rad) 
n 

r 

-66.9 

-65.4 

-60.  8 

(per  rad) 

P 

-5.3 

-20.9 

-33.3 

(per  rad) 

q 

0 

0 

0 

(per  rad) 

*'r 

13.7 

15.2 

17.0 

C.  (per  rad) 

-62.7 

-59.7 

-58.7 

C„  (per  rad) 

g 

D 

0 

0 

C  (per  rad) 

0.3 

0.1 

-0.1 

C  (per  rad) 

0.2 

3.5 

6.  3 
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Table  1.  Continued 
c.  M  ==  1.10,  Altitude  =  10,000  ft 


Trim  Values 

Level 

Flight,  1  q 

Turning 
Flight,  5  q 

Turning 
Flight,  10  c 

(X  (deg) 

0.1 

3.5 

6.0 

4>(deg) 

0 

78.5 

84.2 

(deg/sec) 

0 

7.6 

15.4 

(per  rad) 

-335.7 

-343.3 

-367.3 

(per  rad) 

0 

0 

0 

(per  rad) 

0 

0 

0 

(per  rad) 

0 

0 

0 

(per  rad) 

153.0 

167.4 

184.6 

(per  rad) 

0 

0 

0 

(per  rad) 

-64.6 

-65.1 

-65.3 

(per  rad) 

-0.1 

-3.3 

-6.6 

(per  rad) 

0 

0 

0 

(per  rad) 

10.0 

10.2 

10.5 

(per  rad) 

-55.8 

-56.4 

-57.1 

[per  rad) 

0 

0 

0 

(per  rad) 

0.3 

0.3 

0.3 

(per  rad) 

-1.0 

-1.0 

0.6 
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Table  1 .  Continued 
d.  M  =  1.10,  Altitude  =  40,000  ft 

Trim  Values 

Level 

Flight,  1  q 

Turning 
Flight,  5  g 

Turning 
Flight,  10  g 

a (deg) 

2.6 

9.3 

16.2 

0 (deg) 

0 

79.5 

84.2 

• 

ij,  (deg/sec) 

0 

8.  5 

17.4 

C  (per 

q 

rad) 

-339.7 

-406.0 

-425.8 

a 

rad) 

0 

0 

0 

(per 

P 

rad) 

Q 

0 

0 

(per 

r 

rad) 

0 

0 

0 

Cj^  (per 

q 

rad) 

161.2 

203.4 

216.3 

C  , (per 

rad) 

0 

0 

0 

C  (per 
r 

rad) 

-65.0 

-65.0 

-62.7 

C  (per 
''p 

rad) 

-2.2 

-10.  8 

-19.2 

(per 

q 

rad) 

0 

0 

Q 

Cj  (per 
r 

rad) 

10.1 

10.9 

12.4 

Cj  (per 

rad) 

-56.1 

-58.0 

-53.6 

C,  (per 

q 

rad) 

0 

0 

0 

C  (per 
r 

rad) 

0.3 

0.2 

0 

C  (per 
^P 

rad) 

-0.3 

1.5 

3,4 
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Table  1.  Continued 
e.  M  =  3^0,  Altitude  =  10,000  ft 


Trim  Values 
a  (deg) 


i(>(deg) 

iV<deg/sec) 


C  (per  rad) 
m 
q 

C  (per  rad) 
in* 

CL 

C  (per  rad) 
m 

P 

(per  rad) 

Cj^  (per  rad) 

q 

(per  rad) 
(per  rad) 

C  (per  rad) 
n 

P 

C  (per  rad) 

''q 

C.  (per  rad) 
*r 

Cjj  (per  rad) 
P 

C,  (per  rad) 

q 

C  (per  rad) 
”r 

C  (per  rad) 
^P 


Level 

Flight,  1  g 

0.2 

0 

0 

-227.0 

0 

0 

0 

122.0 

0 

-31.7 

0.7 

0 

4.5 

-40.5 

0 

0.2 

-0.7 


Turning 
Flight,  10  g 

1.1 

84.3 

4.9 

-228.1 

0 

0 

0 

122.4 

0 

-31.8 

0.1 

0 

4.6 

-40.6 

0 

0.2 

-0.6 


Turning 
Flight,  25  g 

2.6 

87.  7 

12.2 

-230.4 

Q 

0 

0 

123.5 

0 

-32.0 

-0.9 

0 

4.6 

-40.8 

0 

0.2 

-0.3 


99 


AEDC-TR-80-11 


Table  1.  Concluded 

f.  M  =  3^0,  Altitude 

=  40,000  ft 

Trim  Values 

Level 

Fliqlit,  1  g 

Turning 
Flight,  10  g 

Turning 
Flight,  25 

s(deg) 

.5 

3.6 

B.4 

4i(deg) 

0 

84.3 

87.7 

Ji;  (deg/sec) 

0 

5.5 

13.7 

(per  rad) 

-227.4 

-232.3 

-252.2 

q 

C  (per  rad) 
ni* 

0 

0 

0 

a 

C  (per  rad] 

0 

0 

0 

(per  rad) 

0 

0 

0 

C  (per  rad) 

122.1 

124.7 

131.4 

Ct  (per  rad) 

0 

0 

0 

(per  rad) 

-31.  7 

-32.1 

-32.4 

(per  rad) 

0.6 

-1.6 

-4.9 

C  (per  radl 

0 

0 

0 

C-  (per  rad) 

4.5 

4.7 

4.9 

*’r 

(per  rad) 

-40.5 

-41,0 

-42.0 

P 

C,  (per  rad) 

0 

0 

0 

q 

C  (per  rad) 

0.2 

0.2 

0.3 

C  (per  rad) 

-0.7 

-0.1 

1.0 

P 
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Table  2.  Yaw-To-Turn  Configuration 


Trim  Values 

a.  M  =  1.3 

Altitude  =  10,000  ft, 
Level  Fliqht,  1  q 

Altitude  =  20,000  ft. 
Level  Fliqht,  1  q 

CL 

6.5 

9.3 

0 

0 

• 

0 

0 

q 

-409 . 7 

-398.7 

0 

0 

p 

0 

0 

0 

0 

0 

0 

% 

0 

0 

r 

-389.3 

-309.7 

P 

-4.1 

-5.1 

C 

n 

q 

0 

^  0 

0 

0 

-4.9 

-5.1 

0 

0 

0 

0 

‘"''p 

0.8 

0.8 
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Trim  Values 
01 


Table  2.  Continued 
b.  M  =  3.0,  Altitude  =  10,000  ft 


Level 

Flight,  1  q 

1.  2 

0 

0 

-298.0 

0 

0 

0 

0 

0 

-309.7 

-0.8 

0 

0 

-3.2 

0 

0 

-0.6 


Turning 
Flight,  5  g 

5.3 

78.  5 

2.8 

-287.7 

0 

0 

0 

0 

0 

-307.0 

-0.5 

0 

0 

-3.6 

0 

0 

-0.2 


Turning 
Flight,  10  g 


9.1 


84.3 


5.7 


-302.9 


0 

0 

0 

0 

0 

-244.5 


-1.2 


0 

0 


-3.7 


0 

0 


-0.1 
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Trim  Values 
a 


Table  2.  Concluded 
c.  M  =  3.0,  Altitude  =  40,000  ft 


Level 

Flightr  1  g 


Turning 
Flight,  5  q 


3.9 


13.6 


m* 


in. 


0 

0 

-287.9 

0 


0 

0 

0 

-310.0 

-0.5 

0 


-3.5 


7B.5 

3.2 

-303.0 

0 

0 

0 

0 

0 

-182.7 

-3.0 

0 

0 

-3.7 

0 

0 


-0.2 


0.1 


Turning 
Flight,  10  q 

22.5 

84.3 

6.4 

-301.8 

0 

0 

0 

0 

0 

-172.6 

-6.5 

0 

0 

-4.1 

0 

0 

1.0 
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Table  3.  Mi»ile  Physical  and  Mass  Characteristics 


Bank-to-Turn 

Yaw-to-Turn 

Mass 

5.75  slugs 

43.80  slugs 

Ij^(body  axis} 

0,34  slug-ft^ 

9.4  slug-ft^ 

lylbody  axis) 

34.1  slug-ft^ 

745  slug-ft^ 

I 2 (body  axis) 

34.15  slug-ft^ 

745  slug-ft^ 

Ijjgfbody  axis) 

0.0  slug-ft^ 

0.0  slug-ft^ 

s 

0.458  ft^ 

1.54  ft^ 

d 

0.7633  ft 

1.4  ft 

1 

8.9  ft 

14  ft 

eg 

50  percent  body  length 

50  percent  body  length 
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APPENDIX 

EQUATIONS  DEFINING  THE  TOTAL  AERODYNAMIC  DATA 
IN  THE  STABILITY  AXIS  SYSTEM 

Longitudinal 

Fjj  =  qS [ Cjj(a,  8p,  M)] 

Lateral-Directional 


where 

6  =  Pitch  Control 

P 

=  Roll  Control 
6^  =  Yaw  Control 
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Bank-to-Turn  Missile 


Yaw-to-Turn  Missile 


Data  Matrix  Variables 


a  -6  to  26  deg 

M  0.8  to  3.5 

&  -10  to  20  deg 
P  ® 


a  0  to  24  deg 
M  0.6  to  3.0 
<5p  -20  to  20  deg 


dp  -10  to  20  deg 
-20  to  20  deg 
dy  -20  to  20  deg 


Control  Authority 


d 

p 

-20 

to 

20 

deg 

-20 

to 

20 

deg 

6 

V 

-20 

to 

20 

deg 
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NOMENCLATURE 

(Note:  All  aerodynamic  data  are  referenced 
to  the  stability  axis  system.) 

Drag-force  coefficient,  drag  force/qS 
Lift-force  coefFicient,  lift  force/qS 

Eterivative  of  lift-force  coefficient  with  respect  to  pitch  rate,  dCt/d(qd/2V),  per 
radian 

Derivative  of  lift-force  coefficient  with  respect  to  d,  dCi_/d{a<i/2V},  per 
radian 

Derivative  of  rolling-moment  coefficient  with  respect  to  roll  rate, 
dC(/3(pd/2V),  per  radian 

Derivative  of  rolling-moment  coefficient  with  respect  to  pitch  rate, 
3C(/d(qd/2V),  per  radian 

Derivative  of  rolling-moment  coefficient  with  respect  to  yaw  rate, 
3Cf/d(rd/2V),  per  radian 

Derivative  of  rolling-moment  coefficient  with  respect  to  d,  3C/£l(dd/2V),  per 
radian 

C(Q  Derivative  of  rolling-moment  coefficient  with  respect  to  angle  of  sideslip, 

per  radian 

C(^  Derivative  of  rolling-moment  coefficient  with  respect  to  /3,  dC|/3(/3d/2V), 

per  radian 

CffiR  Derivative  of  rolling-moment  coefficient  with  respect  to  6r,  dCf/^fiR,  per  radian 

Derivative  of  rolling-moment  coefficient  with  respect  to  6y,  per  radian 

Cm  Pitching-moment  coefficient,  pitching  moment/qSb 

Cmp  Derivative  of  pitching-moment  coefficient  with  respect  to  roll  rate, 

3Cn,/3(pd/2V),  per  radian 

Cniq  Derivative  of  pitching-moment  coefficient  with  respect  to  pitch  rate, 

3Cm/3(qd/2V),  per  radian 


107 


AEDC-TR-80-11 


Cn,|,  Derivative  of  pitching-moment  coefficient  with  respect  to  yaw  rate, 

3Cni/3(rd/2V),  per  radian 

Cm^  Derivative  of  pitching-moment  coefficient  with  respect  to  d,  aCm/3(dd/2V), 

per  radian 

Derivative  of  pitching-moment  coefficient  with  respect  to  0,  dCm/dOSd/ZV), 
per  radian 

Cn  Yawing-moment  coefficient,  yawing  moment/qSd  about  missile  eg 

Cnp  Derivative  of  yawing-moment  coefficient  with  respect  to  roll  rate, 

3Cn/fl(pd/2V),  per  radian 

Cnq  Derivative  of  yawing-moment  coefficient  with  respect  to  pitch  rate, 

3Cn/a(qd/2V),  per  radian 

Cor  Derivative  of  yawing-moment  coefficient  with  respect  to  yaw  rate, 

3C„/3(rd/2V),  per  radian 

C„^  Derivative  of  yawing-moment  coefficient  with  respect  to  d,  9Cn/3(dd/2V), 

per  radian 

Cjiig  Derivative  of  yawing-moment  coefficient  with  respect  to  angle  of  sideslip, 

dCn/d0,  per  radian 

C„^  Derivative  of  yawing-moment  coefficient  with  respect  to  0,  dC„/d(0<l/2V),  per 

radian 

Derivative  of  yawing-moment  coefficient  with  respect  to  6r,  9C„/d6R,  per 
radian 

C„^  Derivative  of  yawing-moment  coefficient  with  respect  to  Sy,  SCn/aSy,  per 

radian 

Cyp  Derivative  of  side-force  coefficient  with  respect  to  roll  rate,  dCy/9(pd/2V),  per 

radian 

Cy^  Derivative  of  side-force  coefficient  with  respect  to  yaw  rate,  9Cy/a(rd/2V),  per 

radian 

Cyg  Derivative  of  side-force  coefficient  with  respect  to  angle  of  sideslip,  dCy/d0,  per 

radian 

Cy^  Derivative  of  side-force  coefficient  with  respect  to  6y,  dCy/d&y,  per  radian 
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eg 

D/R 

d 

Fx  . 
Fy 
Fz 
g 

Ixily.lz 


•xz 


j 

M 

Mx 

My 

Mz 

tn 


p.q.r 


q 

s 

s/p 


V 


W„.  w 
z 


a 


AEDC-TR-BO-11 

Center-of-gravity  location,  percent  body  length 
Dutch  roll 

Reference  body  diameter,  ft 
Force  acting  along  X-stability  axis,  lb 
Force  acting  along  Y-stability  axis,  lb 
Force  acting  along  Z-stability  axis,  lb 
Acceleration  of  gravity,  ft/sec^ 

Moments  of  inertia  about  X-,  Y-,  and  Z-body  axes,  respectively,  slug-ft^ 
Product  of  inertia,  slug-ft^ 

Imaginary  number,  -I 
Mach  number 

Moment  acting  about  X-stability  axis,  ft-lb 
Moment  acting  about  Y-stability  axis,  ft-lb 
Moment  acting  about  Z-stability  axis,  ft-lb 
Mass,  slugs 

Components  of  fl  about  X-,  Y-,  and  Z-body  axes,  respectively,  radian/sec 
Dynamic  pressure,  c  W2,  Ib/ft^ 

Body  reference  area,  ft^ 

Short  period 

Total  velocity,  ft/sec 

Undamped  natural  frequency,  radian/sec 

Damping  ratio 

Angle  of  attack,  deg 

Angle  of  side  force,  deg 
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$p  Pitch  control  deflection,  per  radian 

&R  Roll  control  deflection,  per  radian 

5y  Yaw  control  deflection,  per  radian 

Air  density,  slug/ft^ 

a  Real  part  of  complex  variable 

0  Angle  of  roll,  deg 

SUPERSCRIPT 

Derivative  with  respect  to  time 
K  1,000  of  feet 


no 


